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We present a theoretical analysis of recent experiments on current-driven vortex dynamics in 
the Corbino disk geometry. This geometry introduces controlled spatial gradients in the driving 
force and allows the study of the onset of plasticity and tearing in clean vortex lattices. We describe 
plastic slip in terms of the stress-driven unbinding of dislocation pairs, which in turn contribute to the 
relaxation of the shear, yielding a nonlinear response. The steady state density of free dislocations 
induced by the applied stress is calculated as a function of the applied current and temperature. A 
criterion for the onset of plasticity at a radial location r in the disk yields a temperature-dependent 
critical current that is in qualitative agreement with experiments. 



I. INTRODUCTION 

In the mixed state of type-II superconductors the mag- 
netic field is concentrated in an array of flexible flux 
bundles that, much like ordinarjj_matter, can form crys- 
talline, liquid and glassy phases. Era In clean systems the 
vortex solid melts into a flux liquid via a first order phase 
transition!!] If the barriers to vortex line crossing are 
high, a rapidly cooled vortex liquid can bypass the crys- 
tal phase and get trapped in a metastable polymer-like 
glass phase, much like ordinary window glass .u The diver- 
sity of vortex structures is further increased by pinning 
from material disorder, which leads to a variety of novel 
glasses Bh 

Of particular interest is the dynamics of the vortex ar- 
ray in the various phases and in the proximity of a phase 
transition. In the liquid phase the vortex array flows 
yielding a linear resistivity. In the presence of large scale 
spatial inhomogeneities, the liquid flow can be, highly 
nonlocal due to interactions and entanglement BO The 
correlation length controlling the spatial nonlocality of 
the flow grows with the liquid shear viscosity, which be- 
comes large as the liquid freezes. At a continuous liquid- 
glass transition this correlation length diverges with a 
universal critical exponent.^ In general, probing velocity 
correlations in driven vortex arrays yields information on 
vortex dynamics within a given phase, as well as on the 
nature of the phase transitions connecting the various 
phases. 

As for ordinary matter, the shear rigidity of the vor- 
tex array can be probed by forcing the vortices to flow 
in confined geometries obtained by engineering suitable 
artificial pinning structures, as discussed for instance in 
Ref. H Large scale spatial inhomogeneities can also be 
introduced in the flow, even in the absence of pinning, by 
applying a driving force with controlled spatial gradients, 
as done recently by the Argonne P4P using the Corbino 
disk geometry sketched in Fig. l.ErO In the Corbino disk 
a radial driving current, I, yields a spatially inhomoge- 
neous Lorentz force that decreases as ~ 1/r, with r the 



distance from the center of the disk. For small applied 
currents, the local shear stresses are negligible and the 
vortex array moves as a rigid body. Larger currents (or 
even vanishingly small currents in a glassy solid) result 
in a strong spatially inhomogeneous stress, o(r, I). The 
solid "breaks up" in concentric annular regions flowing at 
different velocity and the response is highly nonlinear .1121 
The voltage profiles measured experimentally by placing 
a series of voltage taps in the radial direction (see Fig. 1) 
reflect the different dynamical correlation lengths in the 
fluid, plastic and elastic regimes. In a given experiment, 
shear-induced plastic slip is observed to occur at different 
values of the applied current in different regions of the 
sample. This is a direct consequence of the controlled 
spatial inhomogeneity of the shear stress. 
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FIG. 1. The Corbino disk geometry. The external mag- 
netic field is out of the page. The electrical current flows 
radially from the inner circumference to the outer rim of the 
disk. The radial current density gives rise to a Lorentz force 
which causes the vortices to moves in circular orbits around 
the disk, without crossing the sample edges. 
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Slippage occurs first in the inner region of the disk where 
the Lorentz force is largest, and propagates to the outer 
regions as the current is increased. In this regime of plas- 
tic response the dynamical correlation length controlling 
velocity inhomogeneities can be identified with the sep- 
aration between free dislocations that are continuously 
created and annihilated by the strong local shear stresses. 

In this paper, we present a theoretical analysis of these 
experiments. A preliminary account of part of this work 
was presented in Ref. [l2|. We describe the onset of plastic- 
ity in the driven vortex solid as a nonlinear effect due to 
the stress-induced proliferation of unbound dislocations. 
The static shear stress from the applied current yields a 
Peach-Koehler force that pulls apart neutral, bound dis- 
location pairs present in the lattice in equilibrium. As 
a result, a finite density, n/(er, T), of free dislocations is 
generated in the solid. Such unbound dislocations can 
then contribute to the relaxation of the shear, resulting 
in a highly nonlinear response to the applied current. 
This mechanism for nonlinear shear relaxation has-been 
studied-Jaefore in the context of superfluid. filmstj, 2D 
solidsEJ'Ej and smectic liquid crystal filmalj. In all this 
cases, the response to a spatially inhomogeneous shear 
was considered. Here we adapt this work to the case 
where the external shear is spatially inhomogeneous, but 
still slowly varying on the length scales of interest. 

Although this paper deals specifically with the on- 
set of plasticity in vortex lattices in the Corbino disk 
geometry, the methods and results presented have a 
more general interest. In particular, the Langevin and 
Fokker-Planck descriptions of dislocation dynamics in a 
spatially-inhomogeneous stress field can be used to de- 
scribe nonlinear shear relaxation in other systems and 
geometries. Our work is for instance relevant to the 
defect dynamics in two-dimensional colloidal crystals, 
where individual point defects can-be created and ma- 
nipulated with optical tweezers Although only va- 
cancies and interstitials have been considered so far, 2d 
colloids, where direct real space and time imaging of de- 
fects is possible, may be especially suitable for studying 
in details the dynamics of dislocations in the presence of 
external stresses. 

We begin in Section II with a brief review of the 
Corbino disk geometry and of the qualitatively differ- 
ent response of a vortex liquid and a vortex solid to an 
applied radial current. In the vortex solid, static elastic 
deformations yield a local shear stress a(r, I), which, for 
the case of free boundary conditions, decreases as ~ 1/r 2 . 
This stress can break bound dislocation pairs present in 
the solid in equilibrium, as shown in sections III and IV. 
For a finite external stress, the total interaction energy of 
a neutral dislocation pair exhibits a saddle point with a fi- 
nite barrier at a value x c ~ l/a(r) of the pair separation. 
At zero temperature, the pair unbinds when the stress is 
large enough that x c (a) ~ ao. At finite temperature un- 
binding takes place when the pair escapes via thermal 
activation over the barrier. Using Langevin and Fokker- 
Planck descriptions of the dynamics of bound pairs, we 



evaluate the rate of thermal escape and use this to esti- 
mate the density of free dislocation, n / (a, T) . The onset 
of plasticity is then identified with the proliferation of un- 
bound dislocations. For free boundary conditions (which 
are appropriate for the Argonne experiments), the ap- 
plied stress is largest near the center of the disk and, for 
a fixed current I, dislocations unbind first near the inner 
rim. We can then evaluate the critical current I p \(r,T) 
where the onset of plastic slip occurs as a function of 
temperature T and of the radial distance from the cen- 
ter, with the result 



I P i(r,T) 



Io(T) 
l+A/r 



(1.1) 



with A = 2R$R%ki(R 2 /R 1 )/(R% - R\) a geometri- 
cal factor. The characteristic current scale Io(T) — 
4cc6et/B exp[— bj (Tm — T) 1 ' 2 ] is the current where shear 
induced proliferation of dislocations occurs across the 
entire sample, with Tm the Kosterlitz-Tbauless melting 
temperature and b a numerical constant .Ej As observed 
in experiments, for fixed T, the current I p \ increases away 
from the center of the disk. Conversely, for fixed r, it is 
strongly temperature dependent and it decreases as the 
melting point Tm is approached from below. 



II. THE CORBINO DISK GEOMETRY 

The Corbino disk geometry was used recently by the 
Argonne group to introduce controlled spatial inhomo- 
geneities in driven vortex arrays and study the onset of 
plasticity is, clean samples, near the first order melting 
transition.E2l The same-geometry has been used before to 
minimize edge cffcctscJ and. also to study the Hall effect 
in type-I superconductorsEil. 

The Corbino disk is an annular superconducting slab 
placed in a magnetic field parallel to its axis (chosen here 
as the z direction), as shown in Fig. |. A current / is in- 
jected at the center of the disk and removed at the outer 
boundary, creating a radial current density, 



J(r) = 



2nrt 



(2.1) 



where r is the radial distance from the center and t the 
thickness of the disk. With the magnetic field applied 
along the disk axis, the Lorentz force density on the vor- 
tices is azimuthal, with 



f L (r) = -B J(r)4> 
c 



(2.2) 



and drives the vortices to move in circular orbits around 
the disk, without crossing the sample edges. Here, 
Bq = no4>o is the average magnetic induction in the sam- 
ple and = z x f is the unit vector in the tangential 
direction. In this geometry, the driving force is spatially 
inhomogeneous and stronger near the center of the disk, 
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falling off as ~ 1/r. It naturally yields spatially inho- 
mogeneous distortions of the vortex array that can be 
probed by placing a series of voltage taps in the radial 
direction. As described below, the spatial dependence of 
the voltage profile in this geometry is qualitatively dif- 
ferent in the vortex liquid and solid phases and provides 
a natural way to study the onset of plasticity. 

As described in Ref. [l^, in the liquid state, vortex mo- 
tion is governed by thrJaydrodynamic equation for the 



vortex flow velocity 



7V - 



rjX7 2 ± v 



(2.3) 



where j(T,H) is the friction and rj(T,H) the shear vis- 
cosity of the vo rtex liquid. The flow profile can be found, 
by solving Eq. ( |2.3| ) with suitable boundary conditions .E£l 
The electric field induced by the vortex motion is then 
obtained from 



E = -B x v 

c 



(2.4) 



If the radial width of the disk is much larger than the 
viscous penetration length, 5 = y/rj/j, the effect of the 
boundaries is negligible and the velocity profile simply 
follows the force profile, with v = fi/7, i.e., v ~ 1/r. 
The voltage drop between the n-th and the n + 1-th elec- 
trodes is then given by 



+ 1 



E{r)dr 



B 2 qI 
7C 2 27rf 



In 



+ 1 



Rn 



(2.5) 



where R n is the radial location of the n-th electrode. This 
logarithmic scaling of the voltage characterizes fluid flow, 
with ncgliblc spatial inhomogeneities caused by viscous 
effects at the boundaries. 

In contrast, a vortex lattice with a finite shear modu- 
lus rotates as a rigid body due to the azimuthal applied 
force, provided the shear stresses due to the local force 
gradients are not too strong. This rigid body rotation 
yields a local velocity that scales as ~ r, with 



n,n+l 



E(r)dr 



Tin 



2c 



n+l 



(2.6) 



where u> is the angular velocity of rotation of the vortex 
lattice, given by 



Bal 



7TC7(i?2 



R\) R\ 



R{ 



(2.7) 



with i?i and R2 the inner and outer radii of the disk, re- 
spectively. Above a certain value of the applied current, 
the local shear stresses become strong enough to break 
the lattice bonds, or generate unbound dislocations. In 
this regime, one obtains a plastic response where the lat- 
tice "breaks" into two or more concentric annular sections 
rotating at different angular velocities and slipping past 
each other due to the continuous generation and recom- 
bination of topological defects. Such a plastic response 
is evident in the experiments by the Argonne groupJ13 



The spatial dependence of the Lorentz force makes it 
possible for all of the types of reponse mentioned above to 
coexist at a given value of the applied current. Near the 
inner radius of the disk the shear stress is very large and 
the vortex array flows like a liquid, with a logarithmic 
scaling of the voltage drop with the contact positions. In 
the middle of the disk the motion is plastic, with disloca- 
tions continuously being created and annihilated. Near 
the outer radius the shear stress is very small and the 
vortex array rotates as a rigid body. Evidence for the 
coexistence of these behaviors was obtained in Ref. |l^ by 
measuring the potential drops across successive pairs of 
electrodes. 

In this article, we describe the onset of plasticity in 
the driven vortex solid as a nonlinear effect due to the 
stress-induced proliferation of unbound dislocations. The 
spatially inhomogeneous Lorentz force gives rise to elas- 
tic deformations of the vortex lattice described by the 
solution of the equation 



fd + ftW = 0, 



(2.8) 



where f e ] is the elastic force density, given by, 

/cu = djaij , (2.9) 



with (Tij the stress tensor, 



aij = 2c 66 Uij + (cu - c 66 )5ijUu 



(2.10) 



and Uij — (djUi + diUj)/2 the symmetrized strain tensor. 
Elastic deformations are described in terms of the two- 
dimensional displacement field, u, and cqq and c\\ are 
the shear and compressional moduli of the vortex lattice, 
respectively. Since dense vortex arrays are practically 
incompressible ( cu >> c§q ), Eq. ( (2.8] ) becomes simply 



c 66 V5_u = -fi(r) . 



(2.11) 



The only nonvanishing component of the displacement 
field is in the tangential direction, u(r) = u^r) <f> 
and yields-a stress a r< j, — 1c§§u r § . For free boundary 
conditional at the inner and outer circumference of the 
disk, [c^uj^Rj = [d r u] r —n 2 = , we obtain 



(r) 



r=R 2 



R\ 



(2.12) 



The spatially inhomogeneous stress given by Eq. ( 2.12| ) 
decreases as 1 /r 2 and gives rise to a Peach-Koehler force 
that can unbind dislocations from bound pairs, yielding 
free dislocations that in turn contribute to the relaxation 
of the applied shear. 
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III. DYNAMICS OF NEUTRAL DISLOCATION 
PAIRS 

Both edge and screw dislocations can occur in the vor- 
tex lattice. The geometry and properties of such disloca- 
tion lines has been discussed for instance in Ref. |24|. For 
simplicity here we assume that the disk is sufficiently thin 
that thermally-induced vortex wandering in the direction 
transverse to the applied field is negligible and the vor- 



tices are essentially rigid rods.c3 In this limit, only edge 
dislocations can occur in the Abrikosov lattice. The dis- 
location lines are aligned with the z direction and their 
Burgers vectors lie in the xy plane. The geometry and 
properties of such rigid dislocation lines are the same as 
those of point dislocations in two-dimensional lattices. In 
particular, the bare interaction energy of an isolated pair 
of dislocations of Burgers vectors and b^ 2 \ located 
at i"i and is given by, 

bW-fo-r^bW-tn-ra) ) 
|ri-r 2 | 2 J 



U^(r 1 -r 2 ) = ^ bW-b^ln 



i r l ~ r 2| 

a 



4-7T 



-^[| b W| 2 + |b( 2 )| 2 ], 
a 



where E c ~ cqqo^ is the core energy per unit length of an 
edge dislocation and the coupling constant Kq is given 
by 



K = 4c 66 (cn - c 66 )t/cn w 4c 66 £ 



(3.2) 



The last approximate equality in Eq. ( |3.2| ) holds for in- 
compressible lattices. The (dimensionful) Burgers vector 
b is defined as the jump in the displacement u upon in- 
tegration around a closed contour, 



du = b . 



(3.3) 



We are interested here in the nucleation of free dislo- 
cations from bound pairs. The interaction energy of a 
neutral pair, consisting of two dislocations of opposite 
Burgers vectors, bW = — b^ 2 ' = — b, is 



U B ( P ) 



47T 



In 



an 



cos 2 6> 



2E r t 



(3.4) 



where p — ri — r 2 and 9 is the angle between p and b. 
The quantity 2E c t represents the energy to create a pair 
of straight edge dislocations of length t at a distance do 
relative to a dislocation- free system. It plays the role of a 
chemical potential for dislocation pairs. The strictly two- 
dimensional limit of point vortices (as opposed to the case 
of rigid vortex lines considered here) is recovered by the 
replacement c§§t — ► c§g, where cqq and c§g are the shear 
moduli of a three-dimensional and of a two-dimensional 
lattice, respectively. 

In general, many dislocation pairs will be present in 
the lattice. The interaction of a given pair is then renor- 
malized by a screening plpud of other pairs. Following 
Kosterlitz and Thouless,E3 this effect can be described in 
terms of a scale-dependent dielectric constant, e(p). Ne- 
glecting the angular part of the dislocation interactiaa, 
which is a marginal perturbation at large length scalescj, 
the effective interaction of a neutral pair is given by 



U Q {p) = 



47T 



dp' 
e(p')p' 



2E r t 



(3.5) 



At equilibrium, in the absence of external stresses, the 
probability per unit area To(p) of finding a pair of dislo- 
cation of separation p is given by 



r (p) 



i 



■ exp 



U {p)/k B T 



(3.6) 



and remains very small away from the Kosterlitz- 
Thoulcss melting transition, Tm = K^a^/ (lQirkB)- Fol- 
lowing Ref. |l6|, the density of pairs on scales p is equiva- 
lently defined by the equation, 



dp 



Kq 

e(p)P 



r (p) , 



with 



-T7 K oal 

An = 



(3.7) 



(3.8) 



a dimensionless coupling constant. Equation (3.7) is 
naturally rewritten in terms of a scale dependent stiff- 
ness, K{p) — Ko/e(p), aad a dislocation fugacity, y(p), 
with [y(p)] 2 — p 4 ro(p).EZl The bare fugacity is simply 
yo = y(ao) = cxp [ — Ec t/ksT] . One can then see im- 
mediately that Eq. (3/7) is nothing but the first of the 
Kosterlitz-Thouless recursion relations, usually written 
in terms of the fugacity as 



dy 

dhip 



2- 



K{p) 



y 



(3.9) 



The scale-dependence of the stiffness arises from the po- 
larization of dislocation pairs and is governed by the sec- 
ond KT recursion relation, 



dK 



d\np 



2ttV 



(3.10) 



It is well known that the KT flow equations yield a low- 
temperature phase with finite long-wavelength dielectric 
constant, separated by a high-temperature phase where 
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e(p) — > oo at large scales and unbound-dislocations pro- 
liferate. The melting occurs at K = 4.E3 

To discuss the nonequilibrium nucleation of free dislo- 
cations from bound pairs under the action of an applied 
stress, we need to study the dynamics of neutral disloca- 
tion pairs as they move under the action of their mutual 
interaction and of the Peach-Koehler force due to the 
external stress field. The position r„ of the Burgers vec- 
tor hj^) is assumed to obey a Langevin equation of the 
form£3 



dr vi 
~dt 



O) 



E dU Q (r u - 
8r ■ 



F 



PK(iy) 



(3.11) 



+riT'{t) , 



where Latin indices denote Cartesian components, p 



(") 



is the mobility tensor, and 77, 
with 



(") 



is a random white noise, 



(VP (t)v { f(t>)) = 2k B T^5^6(t - t') . 



(3.12) 



The force F PK ^ is the Peach-Koehler force from the ex- 
ternal stress and it is given by 



Ff^\v v ) = -e ij a 3k (v v )bt ) 



(3.13) 



Notice that in the case of interest here the external stress 
is spatially inhomogcncous and the Peach-Koehler force 
on the n-th Burgers vector is evaluated at the location of 
the Burgers vector. 

As we are interested in the dynamics of a neutral pair, 
it is convenient to introduce the center of mass and rela- 
tive coordinates of the pair as 



P = ri 

R-CM = 



- r 2 , 
ri + r 2 



(3.14) 



The Langevin equations for the center of mass and rela- 
tive coordinates of a single pair are then given by 



dRcMi 
dt 



IHjCjkbi ou (Rcm + p/2) - <T k i (R C m - p/2) + rfi M (t) , 



and 



dpi 
dt 



Pijtjkh 



(?kl (R. 



CM 



p/2) +a fci (R, 



CM 



P/2) 



dUo(p) 

2pij— \-Vt[t) , 

d Pj 



(3.15) 



(3.16) 



where r) = r/ 2 ) - and r] CM = (rj^ + J7 (1) )/2. As- 



suming that all dislocations have the same mobility, i.e., 
/Ay = pij for all n, we obtain 



(Vi(t)Vj(f)) =2k B Tm j 6(t-t') 
k B T 



(v? M (t)vf M (t')) 



PijS(t - 1') 



(3.17) 
(3.18) 



We now assume that both Ri and the width of the disk, 
R2 — Ri, ar e much larger than do, so that the stress field 
can be considered uniform on the scale ag of the initial 
separation between the dislocations of the pair. We then 
expand the Peach-Koehler force about its value at the 
center of mass of the pair. To leading order, we obtain 



dR, 



CM 



v CM (t) > 



dt 



dpi 
dt 



2pi 



dUo(p) 
dp 3 



6jk(Tkl(R-CM)bl 



(3.19) 
r)i(t) . (3.20) 



The center of mass of the pair performs free thermal 
Brownian motion. We will focus below on the dynam- 
ics of the relative coordinate. 

The Langevin equation for the relative coordinate can 
be simplified if we take into account that climb motion 
is much slower than glide, yielding a strong anisotropy 
in the diffusion of dislocations. In general the mobility 
tensor can be written as 



Pij — P-gtidehbj + Pclimb(Sij — bibj) 



(3.21) 



where p g ndc and p c Umb are the mobility associated with 
the motion of a dislocation in its glide plane (defined by 
its Burgers vector and z) and perpendicular to the glide 
plane (climb), respectively, and b — h/a . Since climb 
can only occur with the cceation of vacancies and inter- 
stitials, Pchmb >> Pghdc The motion of our edge 
dislocations is then predominantly unidirectional, along 
the direction of the Burgers vector. It is convenient to 
separate the relative displacement of the pair, p, into its 
components along the glide and climb directions, 



p = hx + (z x h)y 



(3.22) 



where x = b • p and y = (z x b) • p. The Langevin equa- 
tion for the separation along the glide direction is then 
given by 



dx 
~dt 

where 



and 



2p 



glide 



dUo(p) 
dx 



a CT||(R, 



+ V\\(t). (3.23) 



<7||(Rcm) = bi £ij <7jk (Rcm )bk , (3.24) 



(V\\ (t)v\\ (0) = 2k B Tp glidc 5(t - t') . (3.25) 
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The Langevin equation for the pair separation can also 
be written in terms of the total potential energy of the 
neutral pair in the presence of the external stress as 



dx dU(p) . 



dt 



dx 



(3.26) 



where 



U{p) = U {p) + -Pi[e-i 3 b k + eikbj]ajk(HcM)t ■ (3.27) 

For the geometry of interest here the only non- vanishing 
comp onent of the stress tensor is o r $ given in Eq. 



(2.12). The neutral pairs that can be unbound by this 
stress when climb is forbidden are those with Burgers vec- 
tors along the tangential direction of the disk. Denoting 
simply by r the location of the center of mass of the pair 
relative to the center of the disk and assuming that the 
separation ~ ao of the neutral pair is small compared to 
the width of the disk, the interaction energy of the pair 
can be written as 



U{p) w U Q (p) - a xa(r) , (3.28) 
where a(r) = —<J r ^(r) > is obtained from Eq.( [2.12| ). 




FIG. 2. A pair of dislocations with opposite Burgers vec- 
tors. The Burgers vectors are in the tangential direction. The 
sketch shows the orientation of the pair, but it is not to scale. 
Also shown is the local coordinate system with axis along the 
glide and climb directions of the pair used in the text. 



A plot of the interaction energy U(x,y) as a function 
of the separation in the glide direction, x, for finite y is 
shown in Fig. 3. For simplicity we neglect the angular 
part of the interaction in zero shear and simply use Eq. 
(3.5) for Uq. This approximation does not change the 
qualitative behavior of the results obtained below. For 
a finite value of the applied stress, the interaction has a 
saddle point at a location xo(y) on the positive x axis 
(for oy^ < 0), defined by [dU/dx] x — Xo = and given by 
the solution of 



a ap 2 {) = x Q K(p )k B T , 
where pa(y) = \J x^iy) + V 2 ■ For small 7/,E3 



xo(y) ~Xc~ 2— , 



(3.29) 



(3.30) 



where x c is the solution of the equation 
a Q ax c = K(x c )k B T . 



(3.31) 



At low temperatures, where the coupling constant can be 
replaced by its bare value, we obtain 



K ao 



(3.32) 



The neutral pair can unbind by escaping over the barrier. 
This process creates a pair of free dislocations that then 
contributes to relaxing the applied stress. 




FIG. 3. The interaction energy a neutral dislocation pair 
in the presence of an applied shear stress as a function of the 
pair separation along the glide direction, x, for y — 4ao- The 
value used for the applied stress corresponds to x c — 50ao . 
The dashed line shows the interaction energy of the pair for 

(7 = 0. 

At zero temperature the pair can be said to unbind 
when the saddle point Xq occurs at a length scale of order 
ao. In thick samples, like the one used in the experiments 
by the Argonne group, this may indeed be the relevant 
mechanism for unbinding. In thin films on the contrary, 
we will see that thermal activation plays the dominant 
role. 

To obtain a simple estimate of the critical value of the 
applied stress where unbinding will occur at zero or very 
low temperatures, we use the bare value of the pair bind- 
ing energy and consider a pair with 9 — , i.e., y — 0. The 
saddle point is then given by Eq. ( 3.32| ) and unbinding 
can be estimated to occur when x c (r) ~ ao, correspond- 
ing to a critical stress 



„(r) « K /(4*t) 



(3.33) 



This is the simplest "criterion" for the onset of plastic- 
ity. By solving for r, one obtains the critical radius r p i(J) 



G 



where shear induced dislocation unbinding occurs, yield- 
ing slippage of neighboring annular sections of the vortex 
lattice, 



with 



r P i(D = 



In- I 



1/2 



2i??i?o , / R 2 

■ 111 " 



(3.34) 



where Iq = ccQ§t/ Bq is the current where shear induced 
proliferation of dislocations occurs across the entire sam- 
ple. 

If we use typical parameters of the experiment by the 
Argonne group (i?i sw 35/im, R 2 ~ 350/im, t s» 10/im, 
w 4T, and / » 15to^4 in YBCO), we obtain r pl w 
100/iw. The critical radius r p i(J) which marks the onset 
of plasticity increases with current, indicating that, since 
the shear stress is largest near the inner circumference 
of the disk, "plastic flow" occurs first in annular sections 
close to the axis. This behavior is qualitatively consis- 
tent with the experimental observations. An interesting 
result of our simple model is that, at high fields, the cur- 
rent scale Iq , and therefore r p i(J), are independent of the 
field. 



IV. THERMAL UNBINDING 

At finite temperatures, neutral dislocation pairs can 
unbind by going over the barrier via thermal activa- 
tion. This is certainly the most relevant mechanism 
in two dimensions or in thin samples. The theory of 
stress relaxation via nonlinear unbinding of dislocation 
pairs was developed some time ago by Ambegaokar and 
collaboratorsLj for supecfluid films and by Bruinsma, 
Halpcrin, and Zippeliualia for two-dimensional-crystals, 
and recently applied by Franosch and NelsonEj to de- 
scribe nonlinear shear relaxation in smectic liquid crystal 
films. 

Using standard methods, it is convenient to transform 
the Langevin equation for the pair separation p into— a 
Fokker-Planck equation for the density of pairs, r(p),Eil, 
given by, 



d t T = -diJi , 



where 

Ji = 



dUp 
d Pj 



e jk o-ki(r)bi )T + k B Td j T 



(4.1) 



(4.2) 



is the local current of pairs and r the center of mass of 
the pair. 

Again, we separate the relative displacement p in com- 
ponents along the glide and climb directions, as in Eq. 
( |3.22 ), and neglect climb. Introducing a diffusion con- 
stant for glide, D — /igiide/^sT, and letting u = U /ksT, 
the Fokker-Planck equation can be written as, 



d t T = -d x J x 



(4.3) 



J x (x,y) = -2D[(d x u)T + d x T] 
= -2De- u d x {e u T) . 



(4.4) 



In a steady state, Eq. ( [lg| ) reduces to d x J x = and 
the current of dislocation pairs (which change their sepa- 
ration in the glide (x) direction) is only a function of the 
pair separation in the climb (y) direction. It is given by 



J x (y) = 



(4.5) 



where x is an arbitrary position (with x » a , so that 
the shear stress can be considered uniform). For x — * 00 
we expect that the dislocation pair probability density, 
r, vanishes rapidly, yielding 



J, 



(y) dx'e u( - x '' v ) =2D e u ^T(x,y) . (4.6) 

J X 



The integral on the left hand side of Eq. (4.6) is dom- 
inated by the saddle point and the lower limit of inte- 
gration can therefore be extended to —00. The integral 
is evaluated by the method of steepest descent, with the 
result, 



dx'e u{x '' y) 



2vr 



\[d%u{p)}p=p M (y)\ 



a u(pM(y)) 



(4.7) 



To evaluate T(p) on the right hand side of Eq. fi~6| , 
we first notice that in the absence of external stress 
the density V will assume its equilibrium form, F (p) = 
(l/a^e - " " 3 ) . For small stresses and |x| << xm, we as- 
sume that a local equilibrium "barometric" form holds, 
i.e., T(x 2 ,y) = T{x 1 ,y)e-^ X2 ^-< Xl ^ . The absence 
of climb implies that the integrated pair density at any 
y-cross section remains the same before and after the ap- 
plication of stress, 



dx 2 T (x 2l y) = I dx 2 T(x 2 ,y) (4.8) 
dx 2 T{x u y)e~ [u{x ^ y) - u{xi - v)] . 



The main contribution to the integral in the last equality 
comes from the saddle-point of the full potential u(x,y), 
while the integral in the first term is dominated by the 
saddle point x = of the potential in the absence of shear 
stress. After some algebra, we obtain 



cosh(W) 
V k B T ) 



(4.9) 



Finally, the current is given by, 



Uy) ^ 2D V\[9'Mp)} P = P o(y)\ ^ u(po(v)) (4 



~>™>sh(^) 
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The main approximation used in obtaining this result is 
that of rare escapes over a high barrier. This requires 
K » 1 , which is practically always satisfied as K as- 
sumes its smallest value, K — 4, at Tm- We also must 
have xq >> ao, which imposes an upper limit to the 
value of the external stress. 

The dissociation rate of bound dislocation pairs per 
unit area of the lattice is given by 



R 



dt 



dy I dxT(x,y) 



dyJ x {y) ■ 

(4.11) 



To evaluate this integral, we use a saddle-point approxi- 
mation about y = 0, with the result, 



R 



2D 2/ \ ' 
— TV K X c) 



K(x c ) 



K(x c ) 



(4.12) 



Dislocations of opposite Burgers vectors recombine 
when they come within Xo(y) ~ x c of each other. The 
net production rate of free dislocations is then given byt3 



9n f 
dt 



= R- 



(v)x c n 2 , 



(4.13) 



where n/ is the areal density of free disl ocati ons. The 
second term on the right hand side of Eq. (4.13) is the re- 
combination rate, with (v) ss a^aD /ksT the mean glide 
velocity of a free dislocation under the shearing force 
aao . In the stead y state the density of free dislocations 
is given by n f « ^R/{{v)x c ), or 



,Jf(i=)/2 



\K(x c )] 



3/4 



(4.14) 



For comparison, the density of free dislocations in equi- 
librium, in the absence of external stress, is given by 



n°(T) 



1 



-E a t/k B T 



(4.15) 



The applied shear enhances exp onen tially rif(T,a) over 
its equilibrium value. Equation (4.14) is one of the main 
results of this section. 

To find the explicit dependence of n/ on the applied 
stress, we need to consider the relative importance of 
the two length scales entering the problem: the location 
x c of the saddle point and the Kosterlitz-Thouless cor- 
relation length (T) , which measures the proximity to 
the melting transition below Tm- The KT correlation 
length £_ is defined as the length scale above which the 
scale-dependent interaction K(p) can be approximated 
by its large distance limit, K(p >> £_) sa K(p — oo) = 
4[l + a(T)/4] , with a(T) - (T M - T) 1 ' 2 . It is given by 



£_(T) = a e 



(4.16) 



and diverges as T — > T M . In this regime the fugacity 
is given by y(p » » a{T)p- a( - T ^ 2 . If x c » 
the coupling constant and the fugacity in Eq. ( 4.14 ) can 
be replaced by their large scale values and the density of 
free dislocations is given by 



a(T)x, 



-(2+«(T)/2)_ 



2 2(l+a(T)/4) 

(1 + a(T)/4)V4 



which yields 



n f 



2+a(T)/2 



(4.17) 



(4.18) 



The condition x c >> ^_ or 4:ir<j/K « e^ x l a ^ will 
apply for sufficiently small stress a, not too close to Tm- 
As Tm is approached from below, eventually one, obtains 
>> x c . In this regime one can approximates 



1 



K(x c )n4 1 + . 

L 21n(a; c /a ) 

y{x c ) w [ln(x c /a )] 1 . 
The density of free dislocations is then given by 



(4.19) 



1 2 i 
— [ln(x c /a )] ~ o 2 



(4.20) 



For a Corbino disk of finite thickness t the corr£=. 
sponding Kosterlitz-Thouless transition temperatures 2 ] 
Tm = CQota^l 'inks is very large, so that all temperatures 
of interest are well below Tm ■ In this region, ~ ao 



and the coupling constant and the fugacity in Eq. (4.14) 
can be replaced by their bare values. The density of free 
dislocations induced by the external stress is given by 



n f (T,a)*n° f (T)^ 



,K /2 



K 



3/4 



(4.21) 



At a finite temperature below Tm we can define a con- 
dition for "shear-induced melting" as the value of the 
external stress where x c (a) ~ £_(T). Notice that for 
T << Tm, this condit ion r educes to the one used at 
T = to obtain Eq. (3.33). By solving this for a as 
a function of T, we find that near Tm the critical shear 
stress <7 c (r) where unbound dislocations proliferate in the 
vortex lattice is given by 



a M „ ^0 e -»/(T M -T)^ 
4-7T 



(4.22) 



with b > a numerical constant. The temperature de- 
pendence of this simple estimate is consistent with the 
experimental observation by the Argonne group that the 
applied stress (which is determined by the driving cur- 
rent) required for the onset of plasticity decreases with 
increasing temperature (see Fig. 4 of Ref. |l0|). By com- 
bining the various estimates for the critical shear stress 
where unbound dislocations proliferate, i.e. the vortex 
lattice is "shear-melted" , we obtain the schematic phase 
diagram shown in Fig. 4. 
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FIG. 4. A schematic phase diagram in the (a, T) plane. 
The dashed vertical line is the location of the Koster- 
litz-Thouless melting temperature Tm where thermal unbind- 
ing of neutral pairs occurs in the absence of external shear. 
At a finite shear stress, free dislocations proliferate at the 
lower temperature denoted by the thick continuous-dashed 
line. Even at T — a finite shear stress ~ Kq /4tt is sufficient 
to free bound pairs. The dashed portion of the line separat- 
ing the regions where neutral pairs are bound and the vortex 
lattice responds elastically from the region where free disloca- 
tions proliferate and the vortex arrays flows plastically is an 



interpolation (by eye) between the low and high (T 
estimates obtained in the text. 



J- A 



Finally, recalling that in the Corbino disk geometry 
the external str ess h as the spatially inhom ogcnc ous form 
given in Eq. ( [2.12] ), we can solve Eq. ( 4.22 ) for the 
critical radius r p i(J, T) where shear- melting first takes 
place for a given value of applie d current and temper- 
ature. We find that Eq. (|3.34 ) for the critical radius 
at T = generalizes simply to finite temperatures as 
r p i(I,T) is obtained from Eq. (3.34) by the replacement 
J -> Jo (T) = Jo exp [ - b/(T M - T) 1 ^] . 

In the region below Tm and above <7 c (r), the driven 
vortex array can be described as a lattice with a concen- 
tration nf(a(r),T) of unbound dislocations, moving at 
a steady rate along the direction of their Burgers vec- 
tors. Moving dislocations relax a strain. This provides a 
mechanism for nonlinear stress relaxation that naturally 
yields nonlinear IV characteristics in the superconductor. 

Assuming that the annular width of the disk is much 
larger than ao, so that a hydrodynamic description can 
be used, we define a Burgers- vector density, 



B(r,t) =J2b (v) S(r-r u ) 



(4.23) 



which is related to the local strain in the lattice by 

tikdhWi-j = a Bj , (4.24) 

where uiij is the unsymmetrized strain tensor. It is re- 
lated to the displacement field by Wij — diUj. In the 



presence of unbound dislocations, the displace men t field 
is no longer single- valued, as indicated in Eq. (|3.3| ), and 
local deformations of the medium are more conveniently 
described in terms of the local strain toy . Since free dis- 
locations are always created in pairs, the Burgers-vector 
density is conserved, 



n,Li ■ (I 



(4.25) 



where is the current of the i-ih component of Burgers- 
vector density in the j-th direction, namely 



• u)dr uj 
dt 



(4.26) 



The equation for the displacement field is replaced by an 
equation for the local strain, given by 



d t Wi 



1 

no 



(4.27) 



where j = nov (to linear order) is the number current 
density, with v the local velocity of the vortex array. 

In the Corbino disk geometry of interest here, the only 
nonvanishing component of the dislocation current is in 
the azimuthal direction. Combining Eq. ( 4.27| ) (special- 
ized to a steady state) with the results obtained in section 
IV, we estimate, 



(4.28) 



where we have used Eq. (4.10) evaluated at y = to 
obtain the dependence of the dislocation current on the 
shear stress. The dependence on shear rate is highly non- 
linear as x c ~ l/a(r). The field from flux motion is radial 
and its magnitude is 

E(r) = -Bov+{r) . (4.29) 
c 

The net voltage drop across the disk is given by 

AT / FlrX/r - / ,/r / , /,-' i n(r'\\ l+^(*<= > - 



' E(r)dr - dr dr'[a(r')Y 

Ri JRi JRi 



(4.30) 



If x c >> £_(T), we can replace K by its asymptotic 
value. Since cr(r) ~ J, this immediately gives a nonlinear 
dependence AV^ ~ /3+q/8_ rp ne ex p 0nen ^ a j s nonun i_ 

versal, with a - (T M - T) 1 / 2 . 



V. DISCUSSION 



We have presented a theoretical analysis of recent 
transport experiments on vortex lattices in the Corbino 
disk geometry. In the experiments, the onset of plas- 
ticity corresponds to the onset of nonlincarity in the IV 
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characteristics and an unconventional scaling of the 
age with radial distance from the center of the disk.ErtiJ 
The nonlinear voltage-current scaling is found to be 
AV ~ / 3 + a ( T )/ 8 j with a a nonuniversal exponent that 
depends on temperature. 

In our theoretical model, the onset of plasticity is as- 
sociated with the proliferation of free dislocations which 
break away from tightly bound pairs. In general, two 
mechanisms can be responsible for the nucleation of free 
dislocations from bound pairs: an externally applied 
shear stress and thermal fluctuations. In the absence 
of external shear, the thermally induced proliferation of 
free dislocations is simply the Kosterlitz-Thouless melt- 
ing transition of the lattice. A finite shear stress as the 
one imposed on the vortex lattice by the external cur- 
rent in the Corbino experiment can unbind dislocations 
below the KT transition temperature and even at T = 0. 
Using standard methods of stochastic dynamics, we have 
calculated the density of free dislocations as a function of 
applied stress and the temperature. In contrast to eax= 
lier work on nonlineat-. shear relaxation in superfluid,E£l 
solidliao and smecticO films, here the external stress is 
spatially inhomogeneous. As a result, proliferation of 
unbound dislocations occurs for different values of the 
external current at different locations in the disk. 

The mechanism for the onset of plasticity proposed 
here, namely the proliferation of free dislocations, is cer- 
tainly relevant for thin Corbino disks, where the vortex 
lattice is essentially two-dimensional. In order to com- 
pare our results with the experiments we have assumed 
that even in a Corbino disk of finite thickness dislocations 
are essentially rigid over the thickness of the sample — 
in other words screw dislocations are excluded. In this 
case, the problem becomes essentially two-dimensional. 
The energy cost for creating a dislocation pair is, how- 
ever, proportional to the sample thickness. This raises 
considerably the barrier that bound pairs have to over- 
come to unbind via thermal activation. The finite tem- 
perature calculation may therefore only be relevant in 
very thin samples, while in thick disks the physics is cap- 
tured by the simple T = estimate described in Section 
3. Both at zero and finite temperature, our results, how- 
ever, agree qualitatively with the experimental observa- 
tions, indicating that a simple picture of stress-induced 
dislocation unbinding can account for the spatial depen- 
dence of the onset of plasticity. Specifically, we find that 
(1) "shear-induced melting" starts at the inner boundary 
and propagates towards the outer rim of the disk as the 
driving current is increased, (2) the "critical" current for 
the onset of plastic slip of concentric vortex lattice planes 
decreases with increasing temperature, and (3) the calcu- 
lated critical radius where plastic slip occurs for a fixed 
current is of the same order of magnitude as the value 
measured in experiments. 

A direction for future work is the study of the role 
of dislocation climb, which will of course require to cou- 
ple dislocation dynamics to that of vacancies and inter- 
stitials. In addition, it would clearly be very interest- 



ing to consider dislocation dynamics in the presence of 
quenched disorder. 
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